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electrons of spin 1/2 interacting with localized spins S is studied for the 



The model of Kondo chain with M-fold degenerate band of conduction 



case when the electronic band is half filled. It is shown that the spectrum of 
spin excitations in the continuous limit is described by the 0(3) nonlinear 
sigma model with the topological term with 9 = ir(2S — M). For a case 
■ |M — 2S"| =(even) the system is an insulator and single electron excitations 

at low energies are massive spin polarons. Otherwise the density of states 
has a pseudogap and vanishes only at the Fermi level. The relevance of 
^ ' this picture to higher dimensional Kondo insulators is discussed. 

Is 

PACS numbers: 74.65. +n, 75.10. Jm, 75.25,+z 

The problem of co-existence of delocalized and localized electrons in crystals 
remains one of the biggest unsolved problems in the condensed matter physics. 
The only part of this problem which is well understood concerns a situation 
when localized electrons are represented via a single local magnetic moment (the 
Kondo problem). In this case the local moment is screened at low temperatures 
by conduction electrons and the ground state is a singlet. The singlet formation 
is a non-perturbative process and the relevant energy scale (the Kondo tempera- 
ture) is exponentially small in the exchange coupling constant. It is still unclear 
how conduction and localized electrons reconcile with each other when the local 
moments are arranged regularly (Kondo lattice problem). Empirically Kondo lat- 
tices resemble metals with very small Fermi energies of order of several degrees. 
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It is widely believed that conduction and localized electrons in Kondo lattices hy- 
bridize at low temperatures to create a single narrow band. It is not at all clear, 
however, how this hybridization develops. In particular, it is not clear whether 
the localized electrons contribute to the volume of Fermi sea. If the answer is 
positive, a system with one conduction electron and one spin per a unit cell must 
be an insulator. The available experimental data apparently support this point 
of view: all compounds with odd number of conduction electrons per spin are 
insulators 1 . This class of compounds is now known as Kondo insulators. At low 
temperatures they behave as semiconductors with very small gaps of the order 
of several degrees. The marked exception is FeSi where the value of gap is es- 
timated as ~ 700K 2 . The conservative approach to Kondo insulators would be 
to calculate their band structure treating the on-site Coulomb repulsion U as a 
perturbation. The advantage of this approximation is that one gets an insulating 
state already in the zeroth order in U. The disadvantage is that it contradicts 
the principles of perturbation theory which prescibe that the strongest interac- 
tions are taken into account first. It turns out also that the pragmatic sacrifice 
of principles does not lead to a satisfactory description of the experimental data: 
the band theory fails to explain many experimental observations (see Ref. 2 for 
a discussion). 

In this letter I study a one dimensional model of the Kondo lattice at the half 
filling. I show that the insulating state forms not due to a hybridization of con- 
duction electrons with local moments, but as a result of strong antiferromagnetic 
fluctuations. Due to the local doubling of period of the original lattice the con- 
duction electrons become heavier. An interaction of these heavy electrons with 
spin kinks lead to the formation of massive spin polarons. This scenario does not 
require a global antiferromagnetic order, just the contrary - the spin ground state 
remains disordered with a finite correlation length. I suggest that such scenario 
can be generalized for higher dimensions. Kondo insulators in this case are ei- 
ther antiferromagnets (then they have a true gap), or spin fluids with a strongly 
enhanced staggered susceptibility. In the latter case instead of a real gap there 
is a pseudogap - a drop in the density of states on the Fermi level. The recent 
numerical calculations of Yu et. al. 3 also demonstrate a sharp enhancement of 
the staggered susceptibility in one-dimensional Kondo insulators. 

As a model of one-dimensional Kondo insulator I consider the model of Kondo 
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chain at half filling governed by the following Hamiltonian: 

M 1 

H = X] X] i~ 7}( C t+l,a,a C r,a,a + c ^a,a C r-+l,a,a) + J ( c t,a,a® abCr,a,b) S r ] (1) 
r q=1 z 

It describes an M-fold degenerate band of electrons with spin S = 1/2 interacting 
with local spins S. In what follows I shall use the path integral formalism. The 
path integral representation for spins has been discussed in details by many au- 
thors. I would refer a reader to the book of Fradkin 4 . In the path integral spins 
are treated as classical variables S = Srh (fh 2 = 1); the corresponding Euclidean 
action for the model (1) is given by: 

A = 

J dr{J2[iS J du(fh r (u,r)[d u rh r (u,r) x d T fh r (u,r)]) + c* a>a d T c rj0l>a \ - 

H(c*,c;Sm)} (2) 

The first term is the spin Berry phase responsible for the correct quantization 
of local spins. Since the integrand in the Berry phase is a total derivative, the 
integral depends only on the value of fh on the boundary, i.e. on m(u = 0, r) = 
fh(r),m(u = l,r) = (1,0,0). The introduction of the additional variable u is a 
price one has to pay for the fact that the Berry phase cannot be written as a local 
functional of fh r (r). 

I shall follow the semiclassical approach assuming that all fields can be sepa- 
rated into fast and slow components. The fast components then will be integrated 
out and as a result I shall obtain an effective action for the slow ones. This ap- 
proach is self-consistent if the obtained correlation length for spins is much larger 
then the lattice spacing. In ordinary antiferromagnets this requirement is fulfilled 
only for large spins S » 1. As we shall see later, the Kondo chain is semiclassi- 
cal even for 5 = 1/2 provided the exchange integral is small JM < < 1 . I suggest 
the following decomposition of variables: 



fh r = ak(x) + (—l) r n(x)\Jl — a 2 k(x) 2 , 

{kn) = 0, 

c r = i r ^R{x) + (-i) r i) L {x) (3) 

where \k\a « 1 is the fastly varying ferromagnetic component of the local mag- 
netization. Substituting Eqs.( |3|) into Eq.( 0) and keeping only non-oscillatory 
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terms, I get: 

A = J drdxL, 

L = iS(k[n x d T n\) + ^-[r^c^I + JS(an(x))\Jl — a 2 k(x) 2 ]ijjj + 

2nS x (top — term) (4) 

where 

Top — term — — / drdxe^u {n[dun x 9„n]) (5) 

is the topological term first derived by Haldane 5 . As far as the non-electronic 
part of the action is concerned, my derivation repeats the one presented in Ref. 
4. The interaction of electrons with ferromagnetic fluctuations has been omitted; 
it can be shown that at small JM « 1 it gives insignificant corrections. 

The fermionic determinant is calculated later (see Eq.( ^)). Besides of the 
trivial static part it contains the topological term, but with — M instead of 2S. 
This is what one should expect: this change reflects the fact that local spins 
couple with conduction electrons to give the total spin S — M/2. Substituting 
the expression ( |30|) into Eq.( I get: 

L = iS{k[n x d T n\) + — [{d x n? + (d T n) 2 } + 

2n 

— {JSf In -^(£) 2 + n(2S - M) x {top - term) (6) 

71 JS 

Integrating over fast ferromagnetic fluctuations described by k, I get 



A = — / drdx 
2tt 



v- 2 (d T n) 2 + {d x n) 2 \ + tt{2S - M) x {top - term) (7) 

v ~ 2 = 1 + .pmZ(i/js) (8) 



After the rescaling of the coordinates vr — Xq, x — x\ I get the action of the 0(3) 
nonlinear sigma model with the dimensionless coupling constant 

7TV 71 

9= M = Jm 2 + , 2 ^ - (9) 

V ^ J ln(l/J5) 

This constant is small at JM << 1 which justifies the entire semiclassical ap- 
proach. At \M — 2S\ = (even) one can omit the topological term. In this case 
the model Eq.( 0) is the ordinary 0(3) nonlinear sigma model. This model has a 
disordered ground state with the spectral gap 6,7 

A = Jg- 1 exp[-2n/g] (10) 
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and the correlation length £ ~ J a/ A >> a. If \M — 2S\ = (odd) the topological 
term is essential 5 . The model becomes critical and the correlation functions of 
staggered magnetization have a power law decay. The specific heat is linear 
at small temperatures without requiring, however, the single electron density of 
states to be constant at the Fermi level. 

Now I shall evaluate the fermionic determinant 

D[g] = MTr ln[z 7 ^ + (1 + ij 5 )mg/2 + (1 - ij 5 )mg+ /2] (11) 

where g is a matrix from SU(N) group and m is some constant energy scale (in 
the context of the model ( |T[) m = JS and g = (em)). N = 2 in the original 
problem, but it is worth to do the calculation for general N. I shall study the 
expansion of the determinant ( |ll]) in terms of m _1 Vg. The first terms of this 
expansion are independent of m and survive even at m — > oo. I claim that 
the gradient expansion contains a Berry phase. To prove this point I take a 
route which may seem exotic, but I do not know any better way to get the right 
answer. As a preliminary step I consider the chiral Gross-Neveu model with the 
U (M) x SU (iV)-symmetry described by the following action: 

A = 

J d 2 x{ifj a:a 'y IM d l ,ri ata - ^[{fj aia r} b!a )(fj b> pr} a> p) - fe.aTs^aX^r^a^)] (12) 

The Greek indices belong to the group SU(M) and the Latin ones to SU(N). To 
avoid a confusion I emphasise that this model is not equivalent to the original 
model and I consider it only because the effective action for its low energy exci- 
tations are given by the determinant ( [11]). In order to show that I introduce the 
auxilary field Q ab and decouple the interaction term by the Hubbard- Stratonovich 
transformation. Interacting formally over the fermions we obtain the partition 
function for the tensor Q ab : 

Z = J DQ + DQex P (- J d 2 xL), 

L = ^TrQ + Q - 
2c 

MTr ln^fy + (1 + i^)Q/2 + (1 - i l5 )Q + /2] (13) 

As in the standard [/(iV)-invariant Gross-Neveu model, the effective action has a 
saddle point with respect to Q + Q and fluctuations of det Q are massive. This can 
be shown in the standard fashion. Assuming that the saddle point configuration 
of Q is coordinate independent and as such can be chosen as a diagonal real 
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matrix: Q(x) = diag(Xx, ...A r ), I calculate the density of effective action ( [L3D on 
this configuration: 

M\ 2 I \ I 

Ae ff = U\l/2c + ^ln^} (14) 

where A is the ultraviolet cut-off. The saddle point value of Q, being a point of 
minimum of this function, satisfies the following equation: 

A a /c+^ln^ = (15) 

71 A 

which solution is 

A a = Aexp[-7r/Mc] = m (16) 

The performed calculation suggests that for slowly varying fields Q one can sub- 
stitute Q a b in the Tr In in Eq. ( |l^) by mg a b where g is an SU(N) matrix. In 
other words, the effective action for excitations of the model ( |T^) with energies 
«m coinsides with the fermionic determinant ( II) with m given by Eq.( |I~6|). 



On the next step of the derivation I use the fact that due to the identity 
2 T i T 2 = 1/2 — -Pi 2; where P\% is the permutation operator, the chiral Gross-Neveu 
model ( |T2| ) can be rewritten as the model with current- current interaction: 

A = J d 2 x{if]a, a l^d^ a + AcJ^J^}, 

Ji = (fja^l^rja^) (17) 

where r r are matrices - generators of the SU(M) group. The two models differ 
by a term containing a diagonal scattering. This term does not renormalize and 
therefore is not important. Now I apply to the model ( [T7| ) the non-Abelian 
bosonization procedure suggested by Witten 8 (see also the book 9 ). Namely, I 
rewrite its Hamiltonian in the Sugawara form: 

H = Hu(i) + Hsu(N) + Hsu(M), (18) 
H u{1) = 7i J dx[: J R {x)J R {x) : + : J L (x)J L (x) :] (19) 

2tt Gn r 

Hsu(N) = , N + M) E J ^4 Jr{*)Jr{x) ■ + ■ Ji(x)J l L (x) :] (20) 



i=i 



HsU(M) 



G 



E J dx \ N + M ) { '- J *W J *W : + : J l(*)Jl{x) +4c : J&x)j£(x) :] (21) 

where I have introduced the chiral currents satisfying the Kac-Moody algebra; 
for the group SU(N) the corresponding defenition is 

Jr = VR^atabVRMi 
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A = VtaXbVLM (22) 

where t l are the generators of the SU(N) (spin) group and t)r, tjl are right and 
left components of the Dirac spinor 77. J and J A are the U(l) and the SU(M) 
(flavour) currents defined as in Eq.( but with generators of the corresponding 
algebras. Gn = N 2 — 1 and Gm = M 2 — 1 are the total number of generators of 
the su(N) and the su(M) algebras. Currents from different algebras commute. 
Therefore the Hamiltonian ( [Tjp is a sum of three mutually commuting operators. 
Now notice that the interaction term containing right and left SU(M) currents 
do not affect the spectra of SU (M)-singlets. It is well known that the spectrum of 
Hsu(M) is gapful for the given sign of the coupling constant (see, for example, Ref. 
10). From the previous discussion we can conclude that the gapful excitations 
correspond to fluctuations of det Q. Therefore the spectrum below the gap m is 
described by the rest of the Hamiltonian ( [L8|) , in other words by 

Heff = Hu(i) + H S u(N)', (23) 
H u( i) = n J dx[: J R (x)J R (x) : + : J L (x)J L (x) :]; (24) 

r 2n Gn 

Hsu(N) = J dx {N + M) EC JrWJrW ■ + ■ Jl(x)Jl(x) :) (25) 

The Hamiltonian ( |25|) is the Hamiltonian of the Wess-Zumino-Witten model 
on the group SU(N). Its spectrum is the subsector of the free fermionic spec- 
trum generated by the SU(N) current operators. The model is conformally in- 
variant and exactly solvable 11 ' 12 . In order to extract from these results an ex- 
pression for the determinant ( p|) , I rewrite the model ( |23| ) in the Lagrange 
representation 8-12 : 

^(D = -j^r J d 2 x(d^) 2 (26) 

f M Mr 1 
A SU (N) = J d 2 x{—Tr(d tl g + d ll g) + — J q d^ abc Tr(g + d a gg + d b gg + d c g)} (27) 

where g is a matrix from the SU(N) group {g + g = I, det g = 1). The second 
term in the right hand side of Eq.( [27|) is called the Wess-Zumino term. This term 
is topological. Despite of the fact that it is written as an integral including the 
additional dimension, its actual value (modulus 2mN) depends on the boundary 
values of g(x, £ = 0) = g(x) (g(x,£ = 1) = 0). This property follows from the 
fact that the Wess-Zumino term is proportional to the integral of the Jacobian 
of transformation from the three dimensional euclidian coordinates to the group 
coordinates, and so it is a total derivative in disguise. 
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Now we can recognize in Eq.( \27j). Indeed, its first part (except of the Wess- 
Zumino term) represents the first term in the gradient expansion of Tr ln-term 
in Eq.( Indeed, at small momenta \p\ << m this term is equal to 

M 



J ^xTr^Q+d^Q) {2i 



167rm 2 

Now let us write down the Q-field as follows: 

Q(x) = mg(x)e iVri<t,{x \ 

where g belongs to the SU(N) group. Substituting this expression into Eq.( p8|) 
and taking into account that Trg + d ll g = as a trace of an element of the algebra, 
we reproduce Eq.( and the first term of Eq.( |2"T| ). It comes not entirely 
unexpected that the naive gradient expansion misses the important Wess-Zumino 
term. Being a Berry phase this term requires special care and cannot be derived 
from a gradient expansion. In order to avoid possible calculational difficulties 
I have resorted to the exact solution which gives us the following interesting 
expression for the determinant ( [II]): 

Mm 2 

D[g] = A su{N) [g] + _ln(A/m) + 0(m~ 2 ) (29) 

where Asu(N)[g] is given by Eq.( |27]) and the second term represents the static 
part of the determinant. For the particular case g = (cm) I get 

A SU (2) [0ft)) = ^(d^n) 2 + ix{Top - term) (30) 

Let us discuss the excitation spectrum. I shall do it only for M — 1. In this 
case the original model has a combined symmetry SU(2) x SU(2) (the additional 
S'?7(2)-symmetry arises as a particle- hole symmetry at the half- filling) , excitations 
carry two quantum numbers - spin S, and an isotopic spin /. We have established 
that spin excitations, i.e. excitations with 1 = 0, are described by the nonlinear 
sigma model with the topological term( ^). The leading contributions to the low 
energy dynamics come from antiferromagnetic fluctuations which agrees with the 
results of Ref. 3. The corresponding energy scale ( |Il]) is formally resembles 
the expression for the Kondo temperature. m(J) is larger, however, due to the 
presence of a large logarithm. Therefore the RKKY interaction plays a stronger 
role than the Kondo screening - it also agrees with the conclusions of Ref. 3. The 
topological term can be omitted if \M — 2S\ = (even). In particular, it cancels 
for the most physical case S = 1/2, M = 1. The low lying magnetic excitations 
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are in this case massive triplets, as it is for the 0(3) nonlinear sigma model 6,7 . 
This picture is in a qualitative agreement with the strong coupling limit of the 
model ( |I|). Indeed, for J » 1 the ground state of the Kondo chain consists of 
local singlets. Excited states are triplets and are separated by the gap ~ J from 
the ground state. 

In order to describe the fermionic excitations which have S = 1/2,/ = 1/2. 
It follows from ( [|) that the fermionic fields live in a slowly fluctuating field 
JS{&n(x, r)}. For constant n the electrons would have a spectral gap JS. There 
are states in the gap, however, and they are bound states of electrons and solitons 
of the unit vector field n. Such situation is typical for relativistic field theories 
and for the SU(N) chiral Gross-Neveu model it was discussed in Ref. 13. Let us 
consider a slow static configuration n(x) such that n(— oo) is antiparallel to the 
2-axis and n(oo) is parallel to it. Then a straightforward calculation shows that 
there is an electronic bound state on this domain wall with the zero energy. The 
corresponding wave function is equal to 

poo 

*(x)=Texp[- dy(an(y))]*(0) (31) 

JO 

The total energy of this bound state is equal to the energy necessary to create the 
domain wall, which is the gap for the spin excitations. Therefore in the energy 
interval between JS and m single electron excitations are massive spin polarons. 
It is no longer the case if \M — 2S\ = (odd). The nonlinear sigma model becomes 
critical and supposedly belongs to the universality class of the isotropic S = 1/2 
Heisenberg chain. In the critical phase the spin solitons do not have a fixed scale 
and the corresponding bound states can have an arbitrary small energy. It is 
reasonable to suggest that the single particle density of states in this case has a 
pseudogap on the Fermi surface decaying as a power law: p(u) ~ \ui\ a . 

We can generalize this one dimensional picture for higher dimensions. Suppose 
we have a Kondo lattice in three dimensions. Then at half filling we can achieve 
the insulating state by the periodicity doubling via an antiferromagnetic phase 
transition. Suppose, however, that this transition does not occur, but the system 
is very close to it. Then the low lying excitations are spin polarons. Whether they 
have a gap or not depends on the state of the spin system. For the one dimen- 
sional 0(3) nonlinear sigma model the static magnetic susceptibility is strongly 
enchanced at the antiferromagnetic wave vector q = tt (x(7r) ~ m^ 1 ). In the same 
time it is zero at q = 0. The corresponding measurements in Ce^Bi^Pt^ show 
a pronounced drop in x(0) at low temperatures and the neutron measurements 
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show a spectral gap 1 . It indicates that the ground state in Kondo insulators is 
magnetically disordered. I postpone a discussion of three dimensional problem 
to further publications. 

This work was started during my visit to the Department of Physics of 
Chalmers University in Goteborg (Sweden). I am grateful to Prof. S. Ostlund and 
the condensed matter group in Chalmers for the kind hospitality. I am grateful 
to A. Nersesyan and G. Japaridze for inspirational conversations, for P. Coleman 
for the valuable criticism, to Derek Lee for reading the manuscript and to C. Yu 
for sending me her preprint. 
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